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We study a discrete time spatial branching system on Z'^ with 
logistic-type local regulation at each deme depending on a weighted 
average of the population in neighboring demes. We show that the 
system survives for all time with positive probability if the competi- 
P^ I tion term is small enough. For a restricted set of parameter values, 

CLj ' we also obtain uniqueness of the nontrivial equilibrium and complete 

convergence, as well as long-term coexistence in a related two-type 
model. Along the way we classify the equilibria and their domain of 



C^ ■ attraction for the corresponding deterministic coupled map lattice on 



1. Introduction and main results. An important problem from the field 
^ ' of mathematical ecological modeling is to find plausible stochastic models on 

2l ■ the level of individuals for the time evolution of a "population," say, of ani- 

^Q I mals or plants, which live, move — in the case of plants, we think rather of the 

Q ■ dispersal of seeds — and reproduce in a 2-dimensional space, subject to indi- 

vidual random fluctuations. The mathematically simplest class of stochastic 
models one might come up with, namely, branching random walk and its 
t^ . relatives in which individuals do not interact, are not adequate because in 

dimension 2, they virtually never exhibit stable long-time behavior: it is well 
known (see, e.g., [20] and the discussion there) that they will die out locally 
k> , if the branching is (sub-)critical, and grow locally beyond all bounds if it is 

^ '_ supercritical. 

To describe an "old" population, which corresponds mathematically to 
a nontrivial equilibrium situation, one has to introduce some interactions 
among individuals, which is of course also natural from the modeling per- 
spective. A very drastic solution that is frequently used in so-called stepping 
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2 M. BIRKNER AND A. DEPPERSCHMIDT 

stone models (see, e.g., [19]), in the context of population genetics models, 
is to force the population size, or the population size per deme in a spatially 
extended scenario, to be constant, that is, each birth is exactly matched by 
a death in the population. A number of ecological models have been intro- 
duced and studied rigorously in the context of interacting particle systems 
and probabilistic cellular automata on U^ (see, e.g., [7, 8, 18]). In these mod- 
els the state of the system at discrete or continuous time t is described by 
a function ^^f-.H^ ^ S, where in most cases 5 is some finite set of possible 
types. The interpretation is that a site x is vacant at time t if it{x) = and 
occupied by one individual of type i if ^t{x) = i for some i & S \ {0}. An 
individual at site x changes its type or dies at a certain rate (with certain 
probability in discrete time setting), which depends on the neighborhood of 
X. Reinterpreting the type as occupancy numbers, this class can, in prin- 
ciple, accommodate models with a fixed a priori upper bound on particle 
density. On the other hand, it seems more natural to allow arbitrary popula- 
tion sizes or densities, and introduce a self-regulation mechanism which, for 
example, makes individual reproduction super-critical in presently sparsely 
populated regions and subcritical in crowded areas — accounting for stress 
or competition for resources. Such models with explicit space have been 
studied in the ecological literature (see, e.g., [2, 13]), mostly using computer 
simulations and heuristic arguments (see, e.g., [3, 17] for a comparison of 
different approaches). Recently, some variants of models of locally regulated 
populations have been studied in the mathematics literature and the possi- 
bility of long-time survival in certain parts of the parameter space has been 
rigorously proved for a continuous mass model [1, 10, 11]. 

We add to this literature a variant where particles live in discrete demes 
(arranged on Z ) in nonoverlapping generations, which looks as follows: In 
the absence of competition, an individual has on average m > 1 offspring. 
Due to competition, for example, for local resources, the average reproduc- 
tive success of an individual at position x is reduced by an amount of Xxy > 
by each individual at position y. Here Xxy is a finite range kernel on Z . Thus, 
an individual at x in generation n will have a random number of offspring 
with mean given by 

(1) r"- Zl ^xy^niy)) , 



where ^n{y) denotes the number of individuals at spatial position y in gen- 
eration n, and for r G M, we write r"*" = max{0, r}. In particular, if the occu- 
pancy of neighboring sites is so high that the term in brackets is negative, 
no offspring are generated at site x in this generation. For definiteness and 
simplicity, we assume that the actual number of offspring, given the present 
configuration, is Poisson-distributed with the above mean, and independent 
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for different individuals. Once created, offspring take an independent ran- 
dom walk step according to a kernel p from the location of their mother. 
In this way, our model incorporates individual-based random fluctuations in 
the number and spatial dispersal of offspring. 

A formal specification of the model is given as follows: We assume that 
the motion/dispersal kernel p = {pxy)x,y£Z'' ^^'^ the competition kernel A = 
{^xy)x,y&Z'i satisfy the following conditions: 

(Al) The kernel ipxy)x,yeZ'' — iPy-x)x,yeZ'i is a zero mean aperiodic stochas- 
tic kernel with finite range Rp>l, that is, for all x, y G Z : pxy = for 

ll-*^ y||oo ^ -Kp. 

(A2) 0<Xxy = Ao,y-x, Ao := Aqo > and X^y = for \\y - x\\oo > Rx, where 
For a configuration t] € Mt. and x G Z , define 



(2) 



and 



/(x; 7]) := r]{x) m - Xor]{x) - ^ Xxzr]{z) 
\ z^x ) 



(3) F{x\ii) := Y^ f{y;v)Pyx, 

yeZ'' 

that is, f{y;r]) is the expected number of offspring generated at site y and, 
thus, F(x;r]) is the expected number of individuals at x in the daughter 
generation if the present configuration is rj. Let A^(^'"-), [x,n) E Z x Z_|_ be 
independent standard Poisson processes on M_|_. Given ^„, the configuration 
of the nth generation, S,n+i, arises as 

(4) e„+i(x) = iV(-'")(F(x;Cn)), xeZ^. 

By well-known properties of the Poisson distribution, this definition is con- 
sistent with the intuitive description given above. Note that, technically, this 
model is a "probabilistic cellular automaton" (see, e.g., [5]) with countably 
infinitely many possible states at each site. 

As for all 7] € M^ , we have f{x]r]) < mri{x), for m < 1, one can easily 
construct a coupling of {S,n) with a subcritical branching random walk. In 
that case {^n) becomes extinct in finite time with probability 1 starting 
from any finite initial condition. Our first result states in the case m £ (1,4) 
that if the competition is weak enough, the population, starting from any 
nontrivial initial condition, will survive for all time with positive probability. 

Theorem 1. For each m £ (1,4) andp satisfying (Al), there are choices 
of positive numbers Aq = XQ{m,p) and k* = K*{m,p) such that if Xq < Aq and 
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X^XT^o^Ox ^ f^*^0! then the population survives with positive probability, that 
is, 

Pgo [Vn G N, 3x G Z'^ : Cn{x) > 0] > 

for all ^0 with f{x;^Q) > for some x E Z . Furthermore, conditioned on 
nonextinction 

1 ^ 

n=l 

in particular, the origin (and, in fact, any site x G Z'^) will be occupied at 
arbitrarily large times. 

Note that this result as well as Theorem 3 and Corollary 4 below work 
in any dimension d>l (with threshold values \q,k* depending on d), in 
particular, it establishes the possibility of long-term survival in d = 2. 

The small competition coefficients mean that the system will typically 
be able to maintain a high number of particles per site. In this sense, our 
result concerns a "high density regime." Technically, we follow the natural 
path of a block construction in conjunction with comparison with oriented 
percolation, that might be paraphrased as "life plus good randomness leads 
to more life, so show that bad randomness has small probability." We call 
a space-time point occupied if there are enough particles there and not too 
many in the neighborhood (see Definition 6 for details). The definition is 
such that in the corresponding deterministic model (which is a "coupled 
map lattice" in dynamical systems jargon; see, e.g., [4] for a recent survey 
of this field) 

(5) Cn+i{x)=F{x;Cn), xGZ^n = 0,l,..., 

in which the Poisson variables are replaced by their means, an occupied 
site would after finitely many steps "colonize" its neighbors, that is, make 
them occupied as well. Then we control the probability that this remains the 
case under stochastic perturbation. Choosing small competition coefficients, 
we increase the "typical number of particles" per site in the deterministic 
model. Then we use the fact that the relative deviation of a Poisson random 
variable from its mean is typically small if the parameter is large. Finally, the 
finite range of competition and motion kernels allows to compare the set of 
occupied space-time sites with finite-range dependent oriented percolation 
on a suitable sub- grid of the space-time lattice. 

The method can be adapted to a situation of two competing species to 
show that if in addition to the conditions of Theorem 1 the interspecific 
competition is weak enough, then long term coexistence is possible (see 
Theorem 8). 
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Fig. 1. Function (j) for different values of m. 



The logistic map (t>{x) = x{m — \x) 
deterministic dynamical system 



and especially the one dimensional 



(6) 



Xn+l = (l>{Xn) 



play an important role throughout the paper (see Figure 1 for a sketch of (j) 
for various values of m). For example, in Theorem 1 the restriction to tti < 4 
comes from the fact that otherwise the function (j) would not map the set 
{x G M : (j){x) > 0} into itself. The function (f) has two fixed points, namely, 
and {m — 1)/A. For m £ (1,3), it is well known that is repelling and 
{m — 1)/A is attracting, that is, if xi / 0, then the sequence (xn) converges 
to (?Ti — 1)/A, whereas for m>3, there are no stable fixed points. It is in- 
teresting on its own that the former fact can be generalized to the coupled 
map lattice (5) which is a spatially extended version of (6): If the competi- 
tion coefficients are small enough, the fixed point rj = {m — 1)/ ^^ Aqx of the 
function F is globally attracting in the sense of locally uniform convergence. 

Theorem 2. Let me (1,3), p, A satisfying (Al) and (A2) be given. 
Then there exists a positive number k* = K*{m,p) such that if J2x^o "^Ox < 
K*Xo and f{x;Co) > for some x E Z , then (^.„) converges locally (i.e., 
pointwise w.r.t. ze'L'^) to ("i — 1)/X]x'^0a:- 

Notice that, under the assumptions of Theorem 2, we obtain a complete 
classification of the equilibria of (5) and their domains of attraction: If (^„,) 
does not hit the all zero configuration G Z^ after the first step, it is 
attracted by r/ = (m — 1)/X!a; -^Ox- 

Obviously G Z^ is an absorbing state for {^n)-, so the Dirac measure 
in this state is an invariant distribution for (fn). In view of Theorem 1, it 
is natural to ask if there exist nontrivial stationary distributions, and one 
might expect that if the process does not go extinct, its distribution con- 
verges to some unique invariant distribution. A powerful method to address 
this problem is coupling. Let (^„ ) and {Qi ) be versions of the process {^n) 
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introduced in (4). Let N^'''''\ Tvf '"^ and N^'''''\ {x,n) G Z"' x Z+, be inde- 
pendent standard Poisson processes. We define the coupling of {^n ) and 
{^n ) as follows: 

+ ivf."+i)((i.(^;^(2))_i.(^;^(i)))+). 



(7) 



Theorem 3. Let me (1,3), andp, A as in (Al), (A2) be given. There 
are Ag* = Ao*(m,p) > and k** = K**{m,p) > such that if Aq < Aq* and 
J2x^0^0x ^ ^**'^0; then, conditioned on nonextinction of both populations, 

the coupling of {^n ) and {^n ) is successful in the sense that for each finite 
A C Z , there is a random time T, such that 

^W (-r) = ^^2) (^^) fg^ allxeA and n>T. 

Obviously we have Aq* < Aq, k** < k* . We do not know if in the case 
m G (1,3) the inequalities are strict (but certainly the bounds obtained in 
the proof of Theorem 3 are much smaller than those obtained in the proof 
of Theorem 1). The proof of Theorem 3 again uses the block construction 
argument together with Lemma 14 which essentially follows from [6] where 
it was used in a similar spirit to prove uniqueness of the nontrivial invariant 
distribution for multicolor systems. 

Corollary 4. Under the conditions of Theorem 3 the process (^„) has 
two extremal invariant distributions. These distributions are translation in- 
variant. Conditioned on nonextinction, (^„) converges in distribution in the 
vague topology to a random measure distributed according to the nontrivial 
extremal invariant distribution, that is, we have complete convergence. 

Remark 5. (i) To our knowledge, we present here the first rigorous 
result showing the possibility of long-time survival in a locally regulated 
population in d = 2 for a particle-based model allowing multiple occupancy 
(but for particular cases in a continuous-time version, cf. [11], Proposition 
6.4, where the competition acts strictly within-deme, and Proposition 7.9, 
where competition and dispersal kernel must be identical). Notice that be- 
cause of the nonlocality of the competition and the discreteness of time the 
system considered here is not monotone (see [14], Chapter II. 2, for back- 
ground). As the mathematically rigorous investigation of spatial stochastic 
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systems with local regulation terms is still in its infancy, we think it is justi- 
fied to study the phenomenon in several mathematical guises. Furthermore, 
many species do live in discrete generations, and it is well known that dis- 
crete time dynamics can have a much richer behavior than their continuous 
time analogues. This shows up in our model as well; see point 4 below. 

Being honest, one has to admit that the results of this paper, as well as 
those in [1, 11] are still too weak to capture many ecologically interesting 
phenomena. Up to now, all the rigorous results are more of a conceptual 
nature, showing that survival respectively coexistence of several types is 
possible if the interaction terms are weak enough, but giving little clues 
about what realistic sizes of threshold values enabling/excluding survival 
or coexistence might be. This stems from the fact that in order to apply 
comparison with finite-range dependent directed percolation, one usually 
has to keep far away from the true critical values. For example, we have 
little rigorous information about properties of the nontrivial equilibrium 
guaranteed by Corollary 4, apart from the fact that its mean is close to 
the deterministic prediction (m — ^)/J2x^0x when the competition terms 
are small. One would suspect that correlations decay exponentially, but we 
have no rigorous proof. 

Thus, the contribution of these mathematical investigations to the ques- 
tion, how a population or several populations arrange themselves in space 
in order to survive in a (ecologically very interesting) situation of scarce 
resources and, hence, appreciable competition is at present rather limited. 
It appears that more powerful mathematical tools need to be invented in 
order to make rigorous progress in this direction. 

(ii) The Poisson offspring distribution in our model is a somewhat artificial 
choice, which helps to streamline calculations, but is not essential for the 
result. To formulate a more general form of the model, one would need 
a one-parameter family of probability distributions (say, indexed by their 
mean) which includes sub- and supercritical distributions. A natural way 
would be to start with a fixed supercritical offspring distribution and then 
superimpose a "thinning" according to the local weighted density. A nice 
feature of the Poisson distribution is that we can in fact think of it in this 
way. Another feature of the Poisson distribution is that the variance of 
the total number of offspring produced at some site x (given the present 
configuration) and its mean are the same. While it is natural for a "branching 
model" to assume that conditional variance and mean of the size of the new 
generation are of the same order, a general class of offspring distributions 
would allow for different proportionality factors. 

(iii) Our results require that Aq, the on-site competition coefficient, is 
(substantially) larger than the total competition with neighboring sites. 
Thus, they apply to a situation where most of the competition is felt by in- 
dividuals within the same "colony." One can think, for example, of colonies 
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arranged on Z and Aq governing a rather strong population regulation inside 
each colony, whereas the competition Aqx, x ^0, with surrounding colonies 
is of a lower order. 

This is certainly a technical condition which is not necessary for survival, 
but which intuitively helps quite a bit because it prevents the occupancy of 
a site from becoming so big that it would "eradicate" its neighborhood in 
the next step. Note that no such condition is necessary for the continuous- 
time continuous-mass result in Theorem 1.5, 2 b) in [10] (on the other hand, 
unlike [10], we do not need the requirement that the range of A must not 
exceed that of p) . 

Simulations suggest that the system may survive also when Aq and Aoxj 
< ||x|| < R\, are the same or similar (but sufficiently small), but occupancy 
numbers will fluctuate much more wildly than in the scenario treated in 
Theorem 1. On the other hand, with a highly asymmetric competition kernel, 
one observes in simulations the appearance of "fronts" of occupied sites 
moving in the direction of smaller A. This might indicate local extinction 
despite global survival when starting from a finite initial population in such 
a case. 

(iv) As the model is in some sense a stochastic version of a spatial system 
of coupled logistic maps, the restrictions on m in our results are inherited 
from the behavior of (6): When m > 4, the one dimensional deterministic 
dynamical system (6) would "live" only on a Cantor-like set, and the tech- 
nique employed in the proof of Theorem 1 would fail. On the other hand, 
simulations suggest that, even in the case m > 4, the random fluctuations 
can "smooth out" the trajectories so that (4) might survive from initial 
conditions which would drive (5) to extinction in finitely many steps. The 
restriction to m £ (1, 3) in Theorem 3 stems of course from the fact that this 
guarantees a unique stable fixed point of the logistic map. It is unclear if 
Corollary 4 would hold in a situation where (6) has periodic orbits: Then, 
one can see in simulations large regions of space which are "oscillating out of 
phase," which might indicate that the system builds up long range structure 
over large time scales. Similar effects have been studied in [9]. 

(v) We note that the "stepping stone version of the Bolker-Pacala model" 
introduced in Definition 1.3 of [10] can be obtained as a scaling limit of a 
sequence of models considered above: Assume that the parameters of the 
A^th model are given by 

mS'^> = 1 + 



N 



, (AT) ^ anXxy 
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where a,M,mxy,Xxy are as in [10], page 191. Let ^q (x) = [Nfi{x)], where 
fi is some finite measure on Z'^, and define X^ (x) := j^^lj^^lix)- Then X^^' 

converges in distribution on Z)[o^oo)(-^/(^'^)) to X, the solution of (5) on 
page 191 of [10], that is, the stepping stone version of the Bolker-Pacala 
model, with 7 = 1. Of course, this is a remark about finite time horizons. 
Deducing results about the steady states of the Bolker-Pacala model from 
our theorems would presumably require a considerable amount of work. 

(vi) Hutzenthaler and Wakolbinger [12] have shown that (at least in 
the case of within-site competition only) the stepping stone version of the 
Bolker-Pacala model from [10] dies out in any dimension if the carrying 
capacity, which would correspond to (m — ^)/J2x'^0x in our model, is too 
small. Similarly, one would expect that our model, even when m E (1,3), 
will die out when Xxy are too large. Simulations suggest that this is indeed 
the case, but we have no rigorous proof. 

The rest of this paper is organized as follows: in Section 2 we provide 
a basic lemma showing how "occupancy" spreads through space and prove 
Theorem 1, in Section 3 we briefly discuss how the results can be transferred 
to a two-species scenario with (weak) interspecific competition. Section 4 
provides results about the deterministic system (5) and proves Theorem 2. 
These results will be required in Section 5, where we prove Theorem 3 and 
Corollary 4. 

To simplify the notation in the proofs, we will use in the sequel a trans- 
formed version of the kernel A, 

(8) Xxy = KJxy, X^y, 

where we assume that J2y^x7xy = 1- That is, we separate the nondiagonal 
part of A into k := J2x^o^Ox, the total "nondiagonal" competition and the 
normalized kernel 70a; = Xqx/i^ {ixx '■= 0). Nevertheless, we prefer to state 
the theorems in terms of Xxy because these have an intuitive biological in- 
terpretation. For rj G M^ , x G Z'^ and k > 0, we write 



/k(x; 77) := r]{x) f m - Xor]{x) - k-Y^ 'JxzVi^ 



(9) 

and 

(10) F{x;r]):=Y f^,{y;r])pyx. 



Note that this is just (2) and (3) in the new parametrization. 
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2. Survival. The value fnix; rf) is the mean number of offspring at site x 
if the present configuration is rj. The maximal (mean) number of offspring 
at one site in one generation will be denoted by 

2 
Tfl 

(11) ml := max /^(O;?]) = — -. 

^mf 4Ao 

If the number of particles at some site x exceeds M\^ := tti/Xq, then, as the 
term in the parenthesis in (2) and (9) is negative, no offspring is produced 
at this site. Furthermore, let us introduce 



(12) 



m(Ao,K) := 



m 



1 



An + K 



and mA(, := m(Ao,0), 



the deterministic equilibrium values when the nondiagonal regulation term 
is K respectively 0. Note that, for rj = m(Ao, k), we have /k(x; rj) = m(Ao, «;) 
and, therefore, rj^x) = F{x;r]) for all x G Z*^. 

Definition 6. Let rj G M^ . For a pair of positive numbers (ei,e2)i we 
will say that a site x is {81^82) -occupied with respect to rj if 

r/(x) E [eimAo,(l-e2)MAj and r]{y) <{l- 82)Mx^,\\x - y\\oc,<R\. 

We will often say that r][x) is (ei,e2)-occupied, or just occupied if there is 
no risk of confusion, meaning that x is (ei,e2)-occupied with respect to 77. 
As advertised earlier, to prove Theorem 1, we compare the process (^„) 
with oriented percolation on a sub-grid of Z x Z+. The main step is to show 
that if a site is (ei,e2)-occupied with respect to some |„, then in a while its 
neighbors will be also (ei,e2)-occupied with high probability. To this end, 
we consider a perturbed coupled map lattice 



(13) 



Cn+l{x) =F{x-Xn)+5n{x), 




-2i- Ms 

Fig. 2. Graph of the function f{z) = z{m — Aqz)"*". 
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where the perturbation 6n is assumed to satisfy 6n{x) > —F{x; C,n)-, such that 
(Cn) is nonnegative. We wih show that under certain additional conditions 
on the perturbation term the system {C,n) has the desired property. Then we 
view the original process {^n) as a perturbed dynamical system and we will 
see that the conditions mentioned above are satisfied with high probability 
if the competition is weak enough. 

Let us now introduce and explain some notation which will be used in the 
sequel. We denote by p^y the re-step transition probability of a random walk 
with kernel p. As mentioned above, our goal is to show that an occupied site 
colonizes its neighbors in a couple of steps and remains itself occupied. In the 
first step the offspring are distributed according to the kernel p. Thus, there 
is in general no reason why an occupied site should remain occupied after 
one step. Let us fix some m G (l,reT,). By the Local Central Limit Theorem, 
the number 

(14) re* = min{j G N:pQ^rh^ > 1 for all x with ||x||oo < 1} 

is finite. We set 

J^ = {{y,j) eZ^xZ+:pi,y>0,0<j< re*} 

C {iy,j):\\y\\oo <jRp,0< j <n*}. 

Suppose that the site is (ei,e2)-occupied with respect to Co and that there 
is no mass at the other sites. Let us also assume for the moment that the 
perturbation term vanishes and that the competition between individuals at 
different sites is zero, that is, k = 0. We set f{z) = z{m — Xqz)~^ , see Figure 2. 
As this function is unimodal, to find the minimum of / on some interval, 
it suffices to consider the values at the endpoints. If for some positive a we 
have z G [aei?7iAo) (1 — ^2)Mx^^], then 

( f{e2Mxo) = e2Mxom{l - £2), aeifhxo > e2Mxo, 

~ I f{aeimxo) = aeirhxomil - aei -\ j, aeirhx^ <e2Mxo. 

This means that the number of offspring at site is at least meirfixg if £1 
is sufficiently small. Then the offspring are distributed in the neighborhood 
according to the kernel p. In this neighborhood the mass is again multiplied 
by at least m (unless the local mass happens to be already large enough) 
and then distributed according to p. Hence, after k steps the mass at a 
site y is larger than or equal to p^yin^eiifixo- The living space of the whole 
population at this time is the A:th timeslice of J^ which is contained in the 
ball of radius kRp. By the definition of re* , after re* steps the mass in and in 
points with norm one reaches or maybe exceeds the level EifhxQ- Thus, these 
sites are occupied at that time if the mass there and in the i?A-neighborhood 
does not exceed (1 — e2)MxQ- 
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We need some additional conditions on the perturbation term. We set 

(16) X = {(y, n) G Z'^ X Z+ : ?i < n*, ||y||oo < n{Rp + Rx)}. 

Consider the assumptions: 

(Bl),, for all (y,n) G X: F{y- C„) + 8n{y) < (1 - e2)Mx,, 
(B2)5,K for all (y,n) G X: F{y;Cn) > K implies |<5„(y)| < 5F{y;Cn). 

Lemma 7. Assume that m and p are as in Theorem 1. For each K > 
and 5 satisfying m{l — S) > rhy 1, there are choices of positive numbers ei, 
£2, Aq and K* such that whenever 

(17) Ao < Aq and k < k*Ao, 
the following holds: 

Co(0) is {ei,e2)-occupied, (01)^2, (B2)^^ are satisfied 

=^ Cn*ix) are {ei, £2) -occupied for all x with \\x\\ao < 1- 

Proof. Let K > be given. We choose £2 > such that 

(18) m{l-5)[l-£2 -]>m and ml<{l-2£2)Mx^. 

V 771—1/ " 

For the second inequality, we need m <4. Then we choose ei > satisfying 

777 -|- 1 7Tl 

(19) Po>"ei < £2 < £2 7 for all (n,y) G ^. 

" 771 771 — 1 

Note that this choice guarantees 

777 -I- 1 

(20) p'f^ m''£i'mxo < £2 "Iaq < £2^^,, for all (ti, y) G J^. 

'^ 777 

By construction of J^, the number J^min = min{77i"pg : (77, y) G =^} is posi- 
tive. Therefore we may choose A^ such that, for Aq < A^, 

£ifhxo-ymm>K. 

Finally, we choose At* such that, for some a satisfying rh<m — a, 

{l-£2)MxoH*X*o<a. 

Let us first consider the case k = 0. We have to show that 

Cn*{x) G [£imXo,{l-£2)Mxo], Iklloo < 1- 

By (Bl)£2, we have 

(21) Cn+i{x)=F{x;Cn) + Sn{x)<{l-£2)Mx, for ah (x, 77) G A-. 
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This means, in particular, Cn*{x) < (1 — e2)AfAo for ||x||oo < 1- To complete 
the proof for that case, we show by induction on n that 

Cn(y)G[Poj;m"eimAo,(l-e2)MAj, < n < n* and {y,n)eJ^. 

By definition of n* , the assertion of the lemma then follows. For n = 0, the 
claim holds by assumption. If it holds for some n<n* , then, first using (20) 
and (15), then (19) and the first part of (18), we obtain 

(1 - S)f{y;Cn) > (1 - 5)/(pS,^"ei"^Ao) 

> (1 - 5)p[J^m"eimAo • mM - Po^m^ei + — j 

Hence, 

(1 - 6)F{y; (n) = E (^ " ^)f(^-^ (n)Pzy 
z&lA 

> E Poz^''^^ei7nx„Pzy = eimAom"+Voj^\ (y,n) G ^. 

In particular, we have F{yXn) > ^irnxQ'rh'^'^^PQy > K for Aq < Ag. There- 
fore, (B2)^^ applies and from the last display we obtain 

Cn+i{y) > (1 - 6)Fiy;Cn) > eimx,m^+'p^+\ 

This concludes the proof of the induction and proves the lemma in the special 
case K = 0. 

Now let us turn to the case k > 0. Assumption (Bl)^^, (17) and the choice 
of K* imply that 

< K E 7xyCn(y) < K*Ao(l - e2)Mxo < a, 

||x||oo <n{Rx + Rp) - Rx,n<n*, 
where a > satisfies m — a> rh. We obtain 

/(')(x; Cn) := Cn{x){m -a- AoCn(x))+ < /^(x; Cn(x)) 

< Cn(x)(m - AoCn(x))+ =: /(")(x;Cn). 

So we can use the same induction as in the diagonal case. For the lower 
bound estimates, we use /'-'^ and for the upper bound estimates, we use 

We set Co = Co and assume that (Cn) is the solution of (13) with the 
perturbation term 
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Thus, {^n) with S,n = Cn caii be considered as a perturbed coupled map 
lattice. 



Proof of Theorem 1 . Recall the definition of the space-time box in 
(16). For {x,n) G Z"' x Z+, we set 

X{x,n) = {My'^') :{yj) £ {x,n) + X}. 

Consider the events 

Aix,n) = {N^y^^Hml^) < (1 - e2)Mx„{y,j) G {x,n) + X} 

and 

iV(?/j)(t) 



B{x,n) 



sup sup 

{y,j)e{x,n)+X t>K 



1 



<M. 



We say that X{x,n) is good if A{x,n) riB(x,n) holds. First we want to 
show that the probability of a good realization can be made arbitrarily large 
by choosing small Aq. It is of course enough to consider the corresponding 
problem in the space-time point (0, 0). As A{0, 0) implies (Bl)^ and -6(0, 0) 
implies (62)^^ on the event ^(0,0) nB{0,0), Lemma 7 yields 

{Co(0)(ei , e2)-occupied} n (^(0, 0) n B{0, 0)) 

C {Cn*{y)i l|y||oo < 1 (ei,e2)-occupied}. 

By translation invariance, the corresponding statement is also true for all 
{x,n) G Z'^ X Z_|_. Furthermore, we point out that X{x,n) and X{x',n') are 
independent if ||x — x'||oo > 2n[R\ + Rp) or \n — n'\ > n* . 

Let A be the number of points in X and let {N{t))t>o be a standard 
Poisson process. Then we have 

P[^(0,0)] = (1 - a{\Q)f where a(Ao) =P[iV(m^J > {I - e2)My,]. 

According to (18), we have m*^ < (1 — 2e2)M\f^. Thus, for some ci > 0, we 
have 



a(Ao)=P 



N(m 



\)' 



1> 



(1-^2)Ma„ 



ni 



Ao 



< 



N{^V 



m 



l>e2 



Ao 



< exp 



1 



ci4 

' Ao 



Furthermore, by standard large deviation results for Poisson processes, for 
some C2 > and sufficiently large K, we have 



P[i?(0,0)] 



sup 

.t>K 



N{t) 



t 



1 



<5 



1 



sup 

.t>K 



N{t) 



1 



>6 



>(l-exp(-C2<5^i^)) . 
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From the proof of Lemma 7, one can see that making K large corresponds 
to making Aq small. Hence, 

(22) F[{A{0,0) n B{0,0)y] <¥[A{0,Oy] +¥[B{0,Oy] < 9{Xo), 

where ^(Aq) < exp(— c/Aq) for some suitable positive constant c = c{p, m, Rx). 
This implies 



P[X(0,0) is good] > 1 - 0{Xo) = 1 - (1 - ^p{Xo)f, 

where p{Xo) = (1 — ^(Ao)"'^' )^. Since p{Xo) converges to one as Aq goes to 
and the range of dependence is finite, it is clear that the good sites percolate 
if Aq is small enough. For example, one can apply a result by Liggett, Schon- 
mann and Stacey (see [16], Theorem 26) to show that, for fixed n, the dis- 
tribution of the random field t{x(x,n) is good} dominates the product measure 
i^p[\q) = 02d Ber(p(Ao)) on {0,1}^ xz+ _ Comparison of the process 
(]l{x(x,n) is good})xeZ''xn*z+ with independent oriented percolation concludes 
the proof. D 

3. A competing species model. In this section we consider two processes 
(^n ) and {^n ), modeling, for example, two different species or genetic types 
living in the same habitat and competing for similar (or the same) resources. 
In the absence of the other type, each of them is a version of the basic process 
described in the introduction, possibly with different parameters. 

Let (Axy Ji^^ygz'*' ^li ^ {1>2}, be translation invariant nonnegative ker- 
nels on Z with finite range Rx- These kernels will determine the intra- 

respectively interspecific competition: The average reproductive success of 

(ii) 
an i-individual at x is reduced by each j-individual at y by A^^ . The evo- 
lution of (^n ,^n ) may then be described as follows. Similar to the single 
species model, we define 



/2 {x; el^) , e) ) = el^) (x) m^ - E Aif d^) (y) - E ^^^^ (y) - 

\ y y J 

y 

y 
where rrii is the mean number of offspring of a type i individual in the 
absence of competition. If N^^'"" , N2''^, {x,n) G Z x Z+ are independent 
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standard Poisson processes on M+ , then, given {^n ,^n ), the configuration 
of the next generation is given by 

We obtain the following about long-term coexistence if the competition 
terms are weak enough: 

Theorem 8. For given mi, m2 G (1,4), p*-^-* andp^"^' satisfying (Al) and 
range R\, there are positive numbers A|, A|, nl, K2 and 7* such that if the 
conditions: 

(i) < Xf < A5, Ey^. Ag) < aP<, i G {1, 2}; 
(ii) E,Aa'),E,Ai^^) <7*min{A(^^\Af )}; 

are satisfied, then both populations survive with positive probability, provided 

that for some x,y (^Z we have fi{x;^Q ,S,^'^\0)) > and f2{y;^Q ,^0 )> 
0. Furthermore, conditioned on survival of both populations, 

1 ^ 
liminf -— y^ ]l,^{i),^, .(2),_,^„-, >0 a.s., 

n=l 

that is, we have local coexistence. 

To prove this theorem, one can essentially use the same argument as we 
have used in the proof of Lemma 7 to reduce the case k > to the case 

K = 0. 

4. Results for the deterministic system. In this section we will prove 
Theorem 2. For clarity of exposition, we start with the "diagonal case" 
K = 0. Let us consider more generally a coupled map lattice (Cn) on Z'^, 
defined via 

(23) Cn+l(x) = J2 9{Cn{y))Pyx, X E Z^ 

where {pyx)x,y£Zd is a translation invariant stochastic kernel with finite range 
satisfying (Al) and g : [0, G] — > [0, G] is a continuously differentiable function. 
We think of the single site function g as having as a repelling fixed point 
and another stable fixed point a G (0, G] which attracts (0, G], that is, for 
any xq G (0, G], the sequence (x„) defined through x^+i = g{xn) converges 
to a. [Thus, in particular, ^'(0) > 1, g{G) > 0.] Then obviously C = and 
(^ = a are fixed points of (23), and one is strongly inclined to believe that in 
this well-behaved scenario there are no others. We will say that a dynamical 
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system (ry„) on Z converges locally to a € M if for each finite A C Z and 
each e > there exists Nq such that 

\r]n{x) — a\ <£ for ah x G A and n > Nq. 

Having been unable to find the result we need in the literature, we provide 
Lemma 9 below. Assume the following: 

(DSl) For each a > 0, there exist sequences (a„) and (/3„) such that < 
"0 < a, ^0 = G, an] a, (3n [a and 5( [«„,/?„]) C [«„+!, /3n+i]- 

Note that this implies the following: 

(DS2) There exists a S (0, a) with the following property: 

If Co(0) G [a,G], then there is A'^o £ N such that Cno{x) £ [^iG^]) ll^^lloo < 
1. 

A proof that (DSl) =^ (DS2) is basically a reformulation of the proof of 
Lemma 7. Note that (DSl) holds true, for example, if we assume additionally 
that g is concave (see, e.g., the construction given in Lemma 12). We refrain 
from pursuing the most general conditions for (DSl), but observe that this 
together with Lemma 9 already yields a proof of Theorem 2 in the diagonal 
case K = 0. 

Lemma 9. If Co{^) £ (0,G] for some x G Z'^ and (DSl) holds, then (Cn) 
converges locally to a. 

In the following we will call the set J\fk{A) := {x G Z'^ : infy^A \\x — y\\oo ^ 
k} the /j-neighborhood of A. If j4 = {x}, then we write Nk{x) for the k- 
neighborhood of x. 

Proof of Lemma 9. Let A be a finite subset of U^. We may assume 
that A is a ball with respect to the sup norm. Let {an) and (/3„) be sequences 
from (DSl). Given e > 0, we choose uq such that [3n — oin < ^ holds for all 
n > uq. According to (DS2), there exist a G (0,a) and ni G N such that 

Cnix) > a ^ Cn+miy) > a for ah y with ||x-y||oo < 1- 

Since and a are the only fixed points of g, g'{0) > 1 and a G (0, a], we have 
g{a) >a. It follows that if for all y in the i?p-neighborhood of some point x 
we have Co(?/) ^ Oj then 



Ci(2;) = y29iCoiy))Pyx > a. 



We set 

A':=A/}j^(no+ni)(A) and A./ = A/}j^(„„_i)(A), i G {0, . . . ,no}. 
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AS a;' 



Fig. 3. A sketch of the nested boxes and the dynamical system in one dimension after 
two steps. The superscripts I and r indicate the left respectively the right-hand side of 
the box. Initially the population size is in [ao,/3o] in each site of A' . After two steps the 
population size is in [ao,/3o] on A \Ai, is in [ai,f3i] on Ai \ A2 and in [a2,/32] on A2. 



Note that Aj = AArp(Aj+i) and A.„o = A. By (DS2), there is some time point 
712 G N such that Cn2{x) > a for ah x G A'. We claim that C„2+„(x) > a for 
all X £ Aq and all n > 0. Indeed, during the next ni — 1 steps from time 112 
on, the mass in the points of the i?p-neighborhood of Aq remain bounded 
from below by a. According to (DS2) by the time 71-2 + ni each point in the 
1-neighborhood of A' is bounded below by a. Hence, we are, in particular, 
again in the above situation. 

For simplicity of notation, we assume that Coi^) > a for all x £ A' . We 
need to show that Cn{x) G [a„Q, /3„q] for all 2; G A and n > uq. First, we check 
inductively that for n = 0, 1, . . . ,no we have the following (see Figure 3 for 
an illustration): 



(24) 



(i) Cn{x) G [an,(3n] for X G A„, 
(ii) Ck{x) £ [ak, Pk] for xG Afc\Afc+i,A: = 0,l,...,n- 1. 



For n = 0, (i) is true by assumption, and (ii) is void. Assume that (i) and 
(ii) hold true for some n < hq, let k G {0, 1, . . . , n + 1} and rr G A^ \ A^^i, 
respectively x G A„+i if /c = n + 1. As A/c_i = A/'_Rp(A/c), we have 

Cn+l{x) =^g{Cn{y))Pyx = Yl 9{Cn{y))Pyx G [afc,/3fc] 

y yeAfe-i 

by (DSl), proving (i) and (ii) for n+ 1. 
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To conclude the proof, note that, by the argument above, the set of con- 
figurations C such that 

C(x) > a for X G A', 

C(x) G [afc,/3fc] for x G Afc \ A^+i, /c = 0, 1, . . . ,no - 1, 

C(x) G [a„o,/3„J for x G A„o 

is invariant under the dynamics (23), hence, we have, in particular, for n > 
no, 

Cnix) G [a„o,/3„o] for x G A(= A„o). D 

For the "nondiagonal" case k > 0, we need three more lemmas. Note that 
we only need to consider the case Aq = 1- Otherwise, consider ( defined 
by Cnix) = AoCn(a^)i which solves the iteration given by (9) and (10) with 
Aq replaced by 1 and k by k/Aq- Until the end of this section we write 
^i,K = ?Ti(l, k), m* = ml = m^/4 and fh = fhi = m— 1 [see (12)]. 

Lemma 10. There exist positive n* and 5 such that for k < k* exist 
sequences (an), (Pn) in [niifi — S,nii.o + S] satisfying the following: 

1. a„T?Tii,K, Pnirhi^^,; 

2. If C{y) ^ [an, Pn] for all y € Mr^{x) , then f^{x;C) £ [an+i,Pn+i]- 

Proof. For fixed x G Z'^, we may consider the mapping ( i— > /^(x; C) as 
a function of the restriction of ( to the -RA-neighborhood of x [viewed as an 
element of M where k is the number of points in J\fRs^{x)]. We denote by 
mi,re the vector of length k with all entries equal to fhi^K and by Bs{rni^K) 
the (5-neighborhood of mi,K with respect to sup norm. 

The gradient of C "-^ fni^; C) is given by [we assume that the positive part 
appearing in (9) is not 0] 

d({x)fK{x; C) = m- 2C(x) -kY^ IxyCiy), 

yi^x 

dc,{y)h{x; C) = -t<-lxyC,{x) for y^x. 

Choose positive e, 5 and «;* satisfying 

(|m-2| + 25 + K*(5 + m-l))^<l-e and 
(25) 

.^ ■ \ ^ V^ 

K < min< 



,m-l' ^2(2 + 5), 
For C. G Bs{mifi), we have 

d(:{x)fK,{x] C)<m- 2(mi,o - ^) - k^^ lxy{niifl - 5) 

y^x 

= 2-m + 26 - K{m - 1) + kS 
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and 

dc(x)f^{x\ C,)>m- 2(mi,o + 5) - k ^ 7xy(?"i,o + 5) 

= 2 — m — 26 — n{m — 1) — k6, 
hence, 

\dcix)Ux; 01 < |m - 2| + (m - 1)k + 2d + K6 
and due to (25), we obtain, for k< k* 

(26) (%,)A(x;C))'<l-e. 
For y / X, we have 

|5cfe)/«(y; C)l = K'yxyCix) < Kj^y{mo + 6) 

< K7xj/(m - 1) + S'jxyK < i<i-yxy{2 + 6). 

Consequently, 

(27) ^(a^(,)/.(y;C))'<(2 + <^)V^7^,<(2 + 5)V<|, 

y^x Xyty 

where the last inequality holds if (25) is satisfied. 

Altogether, the above implies that for all C £ Bs{rni^f^) and k < k* we 
have 

(28) l|V/«(x;C)||i<l-|. 

Due to the mean value theorem for all C,,Q' G Bs{mi^Q) exists Q E Bs{mifi) 
such that 

|/.(x;C)-/k(:c;C')I = |V/.(x;C)(C-OI 

<||V/,(x;C)||2-||C-C'l|2<c||C-C'l|2, 



where c = yjl — e/2 < 1. Thus, the claim of the lemma follows. We only need 
to note that /^(x; mi^^) = "t-i.k and that \rhifl — mi^^l < 5 if k < 5/{m — 1) 
which holds by (25). D 

Lemma 11. For each 5 > exists k* > such that whenever k< k* and 
fnix; Co) > for some x £ Z*^, the following holds: For each finite A C Z*^, 
t/iere exists N gN such that Cn{x) G [mi,o — 5, mi,o + 5] for all x £ A and all 
n> N. 
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Proof. Recall our assumption Aq = 1, which implies M^g = m. For all 

X S Z*^, C G [0, m]^ and 6 > 0, we have 

K^^Tx^Cl^;) <mK. 

That implies 

/.,KC(x)) := C{x)im -rriK- C(x)) < /,(x; C) < C(2:)(m - C(x)) =: /„(C(x)). 

The nonzero fixed points of f^^i and /« are respectively rhi = m — rriK — 1 
and TTii^o- Furthermore, if rriK, < 5, then m-i^o ~ ""^z < '^• 

According to Lemma 7, there is ni S N and a > with the property 

(29) Cnix)>a^Cn+niiy)>a, ||a;-y||oo < 1. 

Thus, for each finite A' C Z , there is n2 G N such that Cn2{x) ^ ^ ^or all 
xgA'. 

According to Lemma 12 for each 6 > 0, one can choose k* and sequences 
(a^) and (6„) such that for all k,< k* the following holds: 

Qq < o, 

/K,/([an,&n]),/«([an,&n]) C [an+l,bn+l], 

for some no G N: a„, 6„ E [mo — 5, mo + (5] for all n > hq. 
A construction analogous to the proof of Lemma 9 concludes the proof. D 

The following lemma is a deterministic ingredient in our construction 
[see (DSl)], providing a shrinking sequence of intervals which the one-point 
iteration maps into themselves. Having been unable to find a proof in the 
literature, we provide one here. The property in question will hold for a 
concave / with a repelling and another attracting fixed point and does 
not depend on the particular functional form of /. On the other hand, as we 
also need to consider a slightly perturbed version fs (where in our case the 
perturbation is of a particular functional form), we refrain from generality 
and stick to fs,f'- [0,m] -^ [0, m*], 

(30) fs{x)=x{m — 6 — x)^, f {x) = xijn — x)^ , 

where m* = m?' /A = max/ = f{in/2). Recall that m = ?ti— 1, m^ = m — 5—l 
are the (unique) attracting fixed points of / respectively fs (we think of small 
5). 

Lemma 12. Letting m£ (1,3), consider f,fs as defined in (30). For 
each e > 0, one can choose positive 7 and e, a strictly increasing sequence 
(an), and a strictly decreasing sequence {(3n) with the following properties: 
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(A) There exists Nq G N s.t. a„, /?„ G [fh — e,rh + e] for all n > Nq . 

(B) For all n< No and < 5 < 7: /5(K,/3„]),/(K,/3„]) C [^, %^]. 

Furthermore, Oq > can he chosen arbitrarily small and Pq < m can be 
chosen arbitrarily close to m. 

Proof. We wish to construct the sequences {an) and (/?„) in such a 
way that 

(31) a„ < a,„+i <fh^<m< (3n+i < /5„ 
and 

(32) /^([a„,/3„]),/([Qn, /?.„]) C (a„+i,/3n+i) 
for all n. This together with 

(33) fh — e < hm q„ < hm /3„ < m + e 

n^oo n— >oo 

wiU suffice to conchide, as /^(x) < /^(x) < /(x) for < (5 < 7 and (32) implies 
(B) for each finite Nq and sufficiently small e. The construction is slightly 
different depending on whether the slope of / at its attractive fixed point 
m is G (0, 1), = or G (—1,0), thus, we consider the cases m G (1,2), jtt, = 2 
and m G (2, 3) separately. 

Letting m G (1,2), choose 7 G (0,e) s.t. m — 7 G (1,2). Take arbitrary a^ G 
(0,m — 7) and (3q > m/2 s.t. /^(/So) > f'y{aQ). This guarantees f{[ao,(3o]), 
f.y{[ao,(3o]) C [f^{ao),m*]. Define 

an + Man) . „ 

On+l = ^^ , n>0, 

m* + f f{l3n) + (3n .. 

Pi = 5 and /?„+! = , n>l. 

Note that m* < m/2 in the case considered, so the choice of /?i ensures 
(32) for n = and that /, /^ are increasing on [aoi/Si]- As f-^/{x) > x on 
(0,771^) and f!y{fh^) > 0, we have a„ < On+i < f-yian) for n > 1. Thus, Un /" 
m^. Similarly, observing that x > f{x) > fh for x G {fh, m/2), we have /3„ > 
/3„_l_i > f'y{(3n) for n > 1, hence, /?„ \ fn- This proves (31), (32) and (33) in 
this case. 

Let m = 2. In this case /(m/2) = m* , so the values of f{Pn) cannot be 
decreasing, and we modify the construction as follows: Choose < 7 < e. 
Picking oq G (0,7), define 

fs{an) + an 



an+i 



2 



Pn = -^ + J {2 - 7)74 - Man). n = 0,l,... 
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As above, we have On /" rh-y = 2 — 7. Note that /?„ is the larger root of 
f^{x) = f'y{an), and that the solutions of /^(x) = m^ are m^ = 1 — 7 and 1 
in the case m = 2, so that /(/?«) > f'yiPn) > On+i and (3n \ 1- Hence, (31), 
(32) and (33) are satisfied. 

Finally, let m £ (2, 3). Here, as m > m/2, we need to observe that /([a„, /3„]) 
will contain m* as long as a„ < m/2, so /3„ must not decrease too quickly. 
Furthermore, as f'{rh) < 0, once we come close to m, the roles of the lower 
and upper boundary are interchanged in each step. 

Choose 7 > s.t. m — 7 E (2,3) and m^ > f^{m*) > m/2. Pick uq G 
(0, (m - 7)/2). While (q„ + /^(q„))/2 < m/2, we set 

(Xyi + J'y\Oin) 
"n+l = ^ ■ 

Let no be the smallest integer satisfying (a„o + /^(q„,j))/2 > m/2. We set 

"no +/7(ano) . ^ f^ 



«no+l = 2 2 



(f +/.(-*))• 



Now we choose (3o, . . . ,/3„p s.t. m* < /?, < /3i_i < m and /^(/Si) > Oj+i, ^ = 
l,...,nQ. Note that this is possible because f^{m,*) > ?7i/2. Put /3no+i = 
(/?„„+ m*)/2. 

Let us check (31) and (32) for n < uq: as /^(x) > x for x G (0, m^) and 
f^{m*) < fh.y, the sequence (an)ne{o,...,no+i} is strictly increasing. 
(/9n)nG{o,...,no+i} is strictly decreasing by construction. By definition, we have 

f-yioin) > 2a„+i - a„ > an+i- 

Note that while a„ < m/2, that is, n < no, we always have 

f^i[an,(3n]),f{[an,(3n]) C (a„+i,m*] C (a„+i,^„+i). 

For n > no + 1, define 

(34) a„+i = i(/^(/5„,) + a„), /3„+i = i(/5„ + /(«„)). 
In order to verify (31) and (32) for n > 77,0 + 1, consider 

(35) a£ i—,fh^],b £ {fh,m,) satisfying f (a) <b, f^{b) > a. 

Note that then 

a' = l{a + f^{b)) and b' = ^{b + f{a)) 
fulfill 

a' e{a,fh^),b' e{m,b) and f{a')<b',f^{b')>a'. 
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Indeed, by assumption, we have /^(b) > a, so a' > a. On the other hand, 
f'y{b) < fh^ because f^ is decreasing in [?fi^,m] and h > fh^ = /^(m^). As / 
is decreasing in the considered region, we have 

/(a')</(a)<i(6 + /(a)) = 6'. 

Similarly, b' £ {m,b) and fj{b') > a'. 

Obviously, a = ano+i and b = /3no+i satisfy the condition (35), hence, 
(31) and (32) hold true for n> riQ as well. 

By the above construction, a^ /" a € (?TT./2,m^], (3n\ P £ [fh,m*), where 
{a, (3) solves f{a) = (3, f^{P) = a. For 7 = 0, the unique solution would be 
a = P = fh, for 7 sufficiently small, we have (33). D 

Proof of Theorem 2. Let A be a finite bah in Z'^ and e > 0. Let rii G N 
be such that (29) is fulfilled. For the sequences (a^), (/?n) from Lemma 10, 
choose no s.t. Pn — cin<£ for all n > riQ. Define A' and Aq, . . . , A„g through 

A' = M^no+m){Rx+Rp)i^) and Ai=M^no-i){Rx+Rp)W'^ ^ i^'- ■ ■ ^^o}- 
According to Lemma 11, there exists 712 G N such that Cnix) G [rnifi — 
6,mi^Q + 6] for all x G Aq and n>n2- Then, for simplicity of notation, we may 
assume 722 = 0. Now the rest of the proof is a reproduction of the arguments 
from the proof of Lemma 9. D 

5. Coupling. In this section we prove Theorem 3 and Corollary 4. Let 
us first describe the idea behind the successful coupling. Recall in (7) the 
definition of the coupling (C ,C )• Consider three large (but finite) boxes 
Bi C B2 C B3 C TL and assume that i}^' and (P'' agree on B\ with val- 
ues close to ?TiAo) that they are close to fhx^^ but do not necessarily agree 
on i?2, and that on Bj, all sites are occupied in both systems. In view of 
Lemma 7, we expect that the region of sites which are occupied in both 
systems grows. If the competition is not too strong, the random system 
"follows closely" the deterministic one. Thus, in view of Theorem 2, we 
can hope that the region where both systems are close to the deterministic 
equilibrium "frix^^ is growing as well. Finally, there is a chance that Poisson 
variables whose means are close to each other produce the same realization. 
Therefore, there is also hope that the region where both systems are the 
same grows too. 

Thus, for suitably tuned parameters, we expect that, with high probabil- 
ity, the above situation will reproduce itself after some time on larger boxes 
B[g B'2C B!^. As before, this observation lends itself to a comparison with 
finite range dependent percolation on a coarse grained space-time grid. A 
certain subtlety stems from the problem that the coarse graining must be 
chosen depending on Aq in such a way that the dependence range of the 
percolation does not diverge when taking Aq small. 

For k,len, we set Ak = J\fk(R^+R^+i){0) and ^fe,/ = A/'fc(ij^+ijp+i)+KO)- 
Letting X(y,n), {y,n) G Z'^ x Z+ be the event that for some A^ G N, to be 
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chosen later, the following holds: 



d^)(x)=el^)(x)G 



m — I — 6 771—1 + 5 



An An 



=:I{m,6,Xo) 
for all X £ y + Aj\i, 



(36) C^^\x),4p{x)£l{m,6,Xo) for all x Gy + A^n\An, 
d'^(x),d'^(x) G [£imx„ (1 - e2)Mx,] =: J(m, Aq) 

for all X £ y + Ajn \ A^y . 

Our goal is to show that the process tx(y,n) dominates oriented independent 
percolation on a suitable sub-grid of Z x Z+. The main part of the proof is 
carried out in Lemma 13 below. With this lemma one can, for example, use 
the Liggett, Schonmann and Stacey argument as we have done in the proof 
of Theorem 1. 

Let n* be as defined in (14) and note that this number only depends on 
m and on the kernel p. As we will later choose N large, we will be able to 
choose it as a multiple of n* . In the sequel we will assume that N/n* is an 
integer. 

Lemma 13. For m G (1,3), p as in assumption (Al) and e > 0, there 
exist Aq, K* > such that for each Aq < Aq, k< k,*Xq one can choose N such 
that 

(37) F[X{y,n + N) for ally with \\x - y\\oo < N/n*\X{x,n)] > 1-e 
holds for all x £l/ . 

Proof. Let tti £ (1, 3) and e > be given. Due to translation invariance 
and the Markov property, the left-hand side in (37) does not depend on 
(x,n). Thus, it is enough to prove 

(38) ¥[X{y,N) for aU y with \\y\\^< N/n*\X{<d,Q)] > 1 - e. 
Choose positive e, 5 and k* satisfying 

|m-2| +2(5 + K*((^ + m- 1) < 1 -e and 

* • f <^ ^ 

K < nnn< 



These constants also satisfy (25). Thus, the properties of /^ [see (9)], proven 
in Lemma 10, are preserved. Note that, unlike the situation in Lemma 10, 
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we do not set Aq = 1 here. Furthermore, similar to (26), (27) and (28), we 
obtain 

£ 



||V/,(x;C)||i<l 



(40) 



ifC(2/)e 



m 



1-5 m-l + 5 



Aq Aq 

We choose k^ such that, for all k>kQ, we have 



for all y gMr^^x). 



(41) 



\Ak\ V 2j - \Ak,\ 



(l-0=:c(.)<l. 



We will assume that N >ko. We set Xq = X{0, 0) and X^ = Xn,i n Xn,2 n 
A'at^s, where 

-^TV,! = {Cn\x)=Cn\^) e /(m,(^, Ao) for ah x G ^3iv}, 

'^'7V,2 = {&(.x),^''^\x) G /(m,5, Ao) for all x G AgiV \ ^37v}, 

'^Af.s = 1?)^ (2;),^)^ (x) G J(m, Ao) for all x G ^TTV.Af/n* \ ^6Ar}- 
Furthermore, we define for each n< N the event ^„ by 



^n=ly{x,k) G U ^4(7V~,) X {j}:Cl'\x),ek\x)Gl{m,6,Xo) y 

As A'at implies that X{y,N) holds for all y with ||y||oo < iV/n*, P[A'Ar|A'o] 
is a lower bound for the left-hand side of (38). Therefore, it suffices to show 
F[X^\Xo]<£. Because 

¥[X^\Xo]<F[X^^n^N\'Vo]+n'^^2(^'^N\'Vo] 
(42) 

+ ¥[X^^s\Xo]+¥[-^j,\Xo], 

it suffices to estimate each of the summands. To do this, we will repeat- 
edly use large deviation estimates for Poisson random variables. There are 
constants ci and 5i such that 



(43) 



P[*^|^o]<A^|^4iv|exp 



An 



Now let us consider the first term on the right-hand side of (42). We denote 



r(^.0 . 



(1) c(2)^ 



J'm = ai{Np'>:x£ Z^/ < m},eo>,Co')- We have 

-^ Y: E[ieS)(x)-4^)(x)|i^,|.F^_i] 



<i*A.-i^ E n\CN\x)-c^^\x)\\j^N-i 



^(2), 



iSAf 



xeA^M 
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<u, 



i-3N\ 



E J2 py^ H |v./.(2/;e)iieJv-i(^)-4ii(^)i 



,(2) 



2:e Agjv yeAfRp (x) zeAfR^ (y) 



1 



I A 



SAT I 



J2 \^zfK{y;i)\ J2 P^y 



2GA3JV+1 ^^^""^ ye^fR^{z) xeAsjv 



< ]l-t^„ic(e) 



Asat] V 
1 



2 \A 



37V+1 



ZGA3JV+1 
(1) f^\ t(2) 



A' 



3^+l|.GA3^+i 



E iejvii(^)-ejvii(-)i- 



We can iterate the above argument to obtain on Xq 



1 



l^aivl 



Y: E[\C\^\x)-e^\x)\l^JJ^o] 



m, 



xGAsjv 



(44) 



<c{e) 



<c{e) 



N 



\A. 



AN\ 



E \e^\^)-e:^\^)\ 



(2), 



z&Aa 



N 



\A. 



47V I 



E 25mAo 



2GA4jv\Aj\ 



, sAr|A4Af \ AatI \Afox- 

c(e) — Hi \ — 257nAo < c(e)"25mAo- 



I A 



4Ar| 



From this, we obtain 



(45) 



X&AzN 



\N 



<c{eY2mx,6\A%N\. 



j(2). 



Note that on Xq for all |x| < R\ + Rp we have ^„ (x) = ^n (2;) for all n < 
N -I. 

To estimate the second term of the right-hand side of (42), let (a^) and 
(/?„) be sequences from Lemma 12 satisfying ao < ei{m — 1) and /3o > (1 — 
£2)171. Let K* be small enough for Theorem 1 and Theorem 2 to apply. 
Let A'^o be the number from Lemma 12 such that, for all n > Nq, we have 
a„/((l — (5)Ao),/9n/((l + '5)Ao) S /(jn, Ao,(^). Recall that in the formulation 
of Lemma 12 we have chosen Aq = 1, but it holds for general Aq- We assume 
-^0 < N- If for all X G A/'i?;^+_R,p(0) we have ^0(2;) G [a„/Ao,/3„/Ao], where ^ is 
a version of the processes considered, then there exist positive constants C2 
and 62 such that, for all n < Nq, we have 



6(0) 



Oin+l Pn+l 



An 



An 



iV(0'0)(F(0;eo))^ 



On+l Pn 



+1 



An 



An 
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C25V 



< 



exp 



Ao 
because F(0;^o) G [an+i/((l - 5)Ao),/3n+i/((l + S)Xo)]. It follows that 

(46) P[;f^2 n ^M^o] < N\AjN \ ^47v| expf-^ 

The upper bound for the third term on the right-hand side of (42) is 
obtained as follows: 

F[X^^3\Xo]=n^x G AjN,N/n* :d'^(x),4'^(x) ^ J(m, Ao)|;fo] 



< 



(47) 



N 
3A: G <i 1, ..., — - 1 
n* 

3x G A-rN,k \ AeN : ^S* (x) or ^£1 (x) ^ J(m, Ao)|;fo 



2iV|A7jv,7V/n*-l\A7vL,. , 

< — t^(Aoj, 

where ^(Aq) < exp(— C3/A0) for some positive C3 is defined in (22). 

Let A^ be the smallest multiple of n* larger than I/Aq. Using the above 
estimates, one can choose some positive c and r G N such that 

p[;f^|;fo]<exp(^-^)iv^ 

The right-hand side goes to zero as Aq goes to zero. Thus, (38) follows. D 

Before we turn to the proof of Theorem 3, we need a result about oriented 
percolation. Let 6 G (0, 1) be given and let A{x,n), {x,n) G Z*^ x Z+ be i.i.d. 
Bernoulli random variables with parameter 9. For k <n, we say that {x,k) 
is connected to {y,n), this will be denoted by {x,k) — > {y,n), if there is a 
sequence x = xq, . . . ,Xn-k = y such that Hxj — 2;j_i||oo < 1 and A{xi,k-\-i) = 1 
for i = 1,. . . ,n — k. Let Cq = {(x, n) : (0, 0) — > {x, n)} be the cluster of the 
origin. We call a space time-point (y,n) Co-exposed if there exists a sequence 
yn,-.-,yo such that yn = y, ||?/fe - yfc-ilU < 1, and {yk,k) ^ Cq, k = 1,. . . ,n. 

The next lemma follows from [6]. The idea behind the proof is as follows: 
With the "usual" percolation interpretation in mind, let us call a site {x,n) 
wet if there is a backward path (x,n) = (xo,n), (xi,n — 1), ... , (xn,0) with 
||xj — Xi_i|| < 1 consisting only of open sites, that is, A{xi,n — i) = 1, i = 
0, 1, ... ,n — 1. Otherwise, the site will be called dry. Lemma 7 in [6] shows, 
using a contour-counting argument, that if 6 is sufficiently close to 1, the 
dry sites do not percolate. In fact, this lemma even obtains an exponential 
bound on the tail of the size of the cluster of dry sites containing a given site. 
The next ingredient is complete convergence for oriented percolation ([6], 
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Lemma 8) : When 9 is close enough to 1 , there is a fixed c > and a random 
A'^o such that on {|Co| = oo}, {(a;,n) : (x,n) wet and ||a;|| < en} C Co for all 
n > ATq. In words, any wet site inside the "cone" {(x, n) : ||x|| < en, n > Nq} is 
also connected to (0, 0) by an open path. Fix c' € (0, c). Assume that {|Co| = 
oo}, consider (y,n) with ||y|| < c'n and n > 2Nq, say. If {y,n) is Co-exposed, 
there must be a backward path (y,n) = {yQ,n), {yi,n — 1), ... , (yn,0) with 
{yi,n — i) ^ Cq. By the above, at least the initial n{c — c')/2 of these sites 
must be dry [for otherwise, they would be in Cq, as they must satisfy \\yi\\ < 
c{n — i)]. Hence, there must be a cluster of dry sites containing a point in 
{{x,n) : \\x\\ < c'n} of size at least n{c — c')/2. By the exponential bound 
on the cluster size distribution and the Borel-Cantelli lemma, this does not 
occur for n sufficiently large. 

Lemma 14. If 6 is sufficiently close to 1, then there is a positive constant 
c such that, for large enough times n conditional on {\Cq\ = oo}, there are 
no Co-exposed sites in {x G Z'^: ||x||oo < en}. 

Proof of Theorem 3. Recall the definition of the event X{y,n) from 
(36). Theorem 1 implies that, conditioned on nonextinction of (^n ) and 

(2) 

(^n ), with probability one, there exist some finite time A'^o such that the 
event X{0,Nq) holds. Therefore, we may assume a priori that .^(0,0) holds. 
^ We set N = [N/{2n*)], B = {{x,n) € Z'^ x Z+|||x||oo <N,n<N}, L = 
NZ'^ and K = NZ^. Then we have 



yd 



X 



U ((«>'^)+^)- 



(a,i/)€LxK 



Let II • \\l be the norm on L defined by ||a||i = ||a||oo/-^- To prove the 
theorem, it is enough to show that for each x* G Z'^ there is time T, such 
that ^n (x*) = ^n (x*) holds for all n>T. Let us fix an arbitrary x* G Z'^ 
and let a* G L be such that ||a* — x*||oo < N. We define a process {r]^) on 
the coarse-grained lattice L x K hy 

Vo{a) = tx{a,o) and Vi^{x) = lx{a,u~N), ^>0. 

Note that lx[a,u-N) = 1 for i^ > ensures that ,^^ (y) = ,^^ (y) holds for all 
(y, k) G (a, i^ — N) + B, because any backward in time path starting in (y, k) 
will at time v — N he inside a + A^, where ^'^^ and ^"^ are the same on 
the event X{a,u — N). In particular, r],y{a*) = 1 implies S^f. (x*) = S,}. (x*) 
for all k G {v — N, . . . , v}. We aim at showing that, for suitable choice of 
parameters, the process {tju) dominates oriented percolation on L x iC. To 
this end, we need to estimate 

P[7?^+jv(/3) = 1, ||a - /3||l < l|??^(a) = 1], 
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whereas, due to translation invariance, it is enough to consider the cor- 
responding probabihty for (aju) = (0,0). By the construction of (r/,y) and 
Lemma 13 for each positive e, one can choose Aq, k and N such that 



P[r/^(/3) = l,||/3||L<l|r?o(0) = l] 



> 



N 

;f(z,iv),||z||oo< — ;f(o,o) 



>l-e. 



From the proofs of Theorem 1 and Lemma 13, it can be seen that for x with 
W^Woo ^ 1 the event X{Nx, N) is independent of the Poisson processes [which 
generate {^n^) and {^n^)] outside the box {(y, /c) e Z"^ x Z+ : A; < N, ||y||oo < 
(8A^ + 2){Rx + Rp)}. Therefore, (iji,) can be considered as M-dependent 
oriented percolation on LxK, where M = 20n*{Rx + Rp) > {8N + 2){Rx + 
Rp)/N . Note that M does not depend on N and Aq. Thus, the fact that we 
need to make Aq small does not affect the comparison. 

Let 9 be close enough to 1 such that Lemma 14 holds. For e G (0, (1 — 
VO)'^), where A = \{{a,u) £ L x K :i^ e {0,N},\\a\\L < M}\, we have 

p[7?^(/3) = 1, m\L < i|77o(o) = 1] > 1 - (1 - Vef. 

As in the proof of Theorem 1, according to Theorem B26 in [15], (r/jy) domi- 
nates the nearest neighbor oriented percolation build from the product mea- 
sure i^g on {0, 1}^^^ . Thus, we obtain IP[|C,,| = oo] > 0, where C,^ C L x K 
is the cluster of the origin generated by (r/j/). By Lemma 14, conditioned 
on {\Crf\ = oo}, there is a time T such that the points {a* ,1^) £ L x K with 
i'>T are not C,;-exposed. 

We claim that, for each n >T, conditioned on {\Crj\ = oo}, we have 

Cn {x*) = ^n (ic*). If we assume the contrary, then there must be a path 
(x*,n) = (a;„,n), (x„_i,n — 1), . . . , (a;o,0) in Z'^x Z4. such that Hxj+i — Xj||oo < 
R\ + Rp and Q (xj) / Q (xj) for alH G {0, . . . , n — 1}. From this path, we 
discard the points (xj, i) for which i is not a multiple of N, thus obtaining for 
some integer k the path {xkN,kN),{x(^i^_i)j\f,{k — 1)N), . . . ,{xo,0). To this 
path belongs a path (a*, {k + l)N), (a^N^ kN), . . . , (aO)0) in Lx K where for 
j G {1, . . . ,k} we choose Oj^ such that {x(^j_i)n, {j — 1)N) G {q.jNiJN) + B 
and ao such that ||ao — xo||oo <N. The assumption means that rjiN{c(iN) = 
for alH G {0, . . . , /c}. This contradicts the fact that (a*, {k + 1)N) is not C,,- 
exposed. D 

Proof of Corollary 4. The sequence {S,n), seen as a sequence of 
random measures on Z'^, is relatively compact with respect to convergence 
in distribution in the vague topology because the expectation of ^n{x) is 
bounded uniformly by mX . 

It is clear that Dirac measure in G Z^ is invariant. If there were two 
invariant distributions assigning probability to the configuration 0, then 
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Theorem 3 would imply that they coincide on finite subsets of Z and, there- 
fore, they must be equal. 

It remains to prove the existence of a limiting invariant distribution fj, 
satisfying /^(O) = 0. Let the initial distribution /io be the product measure 
on Z'^ such that ^o{x) = N^'^^''\f}i{Xo,K)) for ah x E Z'^. Let /i„ be the distri- 
bution of ^n- Then the Cesaro average ^/NJ2n=of^n converges along some 
subsequence {N^} to some measure p,. This measure is invariant for (,^.„) 
(see, e.g., [14], Proposition 1.1.8). 

To show /u(0) = 0, it is enough to prove that the restriction of {^n) to Z 
survives with probability 1. At time 0, each site is occupied in the sense of 
Definition 6 with probability 

P[iV(0'0)(m(Ao,K)) G [eifhx„,{l-e2)A'h,]], 

where £i and £2 are as in the proof of Lemma 7. In particular, at time 0, 
there are infinitely many occupied sites. Again, by comparison with oriented 
percolation, we have P^^ [^n = for some n] = because supercritical perco- 
lation starting from infinitely many wet sites does not die out (see, e.g., 
Theorem B24 in [15]). □ 
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